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I ' The 5'f7(2)-Higgs model, with a single Higgs field in the fundamental representation and a quartic 

^ , self-interaction, has a Higgs region and a confinement region which are analytically connected in the 

• parameter space of the theory; these regions thus represent a single phase. The effect of multiple Higgs 

0^ I fields on this phase structure is examined via Monte Carlo lattice simulations. For the case of > 2 

■ identical Higgs fields, there is no remaining analytic connection between the Higgs and confinement 
P— I, regions, at least when Lagrangian terms that directly couple different Higgs flavours are omitted. An 

' explanation of this result in terms of enhancement from overlapping phase transitions is explored for 

'— I 2 by introducing an asymmetry in the hopping parameters of the Higgs fields, ft is found that an 

^ • enhancement of the phase transitions can still occur for a moderate (10%) asymmetry in the resulting 

O , hopping parameters. 

[ The phase structure of the basic SU(2)-Higgs model {i.e. an SU{2) gauge theory coupled to a single 

■ scalar field with a quartic self-interaction and a quadratic term) is clearly of direct relevance to the Higgs 
. sector of the Standard Model, where discussions of spontaneous symmetry breaking and the transition 

\^ \ from a symmetric phase to a Higgs phase are paramount. The addition of extra Higgs fields occurs in a 
^—1 ■ variety of extensions to the Standard Model, including the minimal supersymmetric Standard Model (see 
: e.g. Ref. Hj). 

\ In the case of a single Higgs, the S'C/(2)-Higgs model has three parameters in its continuum formulation 

Qv^ ■ (Higgs self-coupling, gauge coupling, and Higgs quadratic term); in the lattice formulation of these models 
<^ ■ the three parameters become the Higgs self-coupling A, the gauge coupling /3, and the hopping parameter 
K. For a single Higgs field in the fundamental representation, there exist theoretical arguments that there 
must be an analytic connection in (/3, k. A) parameter space between the Higgs and confinement regions 
of the theory [2] . This behaviour is manifested in lattice simulations as a phase-transition line that begins 
at (/3 = oo, Koo > 0) and terminates at a point (/? > 0, k > Koo), with a general progression of this 
termination point to smaller values of with decreasing A Od]. The resulting analytic connection 

between the Higgs and confinement regions (located in the corner of parameter space toward larger A, 
smaller /?, and larger k) is consistent with Elitzur's theorem [5] which says that a local gauge symmetry 
cannot break spontaneously. For a recent discussion of the breaking of global subgroups within the SU{2)- 
Higgs model and their connection to the phase diagram, see Ref. [6]. We note in passing that Ref. |6] also 
mentions the S'C/(2)-Higgs model where the Higgs field is in the adjoint representation rather than the 
fundamental; in that case the theoretical arguments of [2] do not apply and the center symmetry can be 
broken spontaneously. 

Many basic properties of the phase structure for a single Higgs field in the fundamental representation 
are well understood, though there is no consensus on a detailed understanding of the nature of the phase 
transition (PT) across the entire parameter space. For small A, the PT is demonstrably of first order 
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[7] but the PT strength weakens with increasing A, comphcating the classification of the PT. The most 
effective approaches used to address this issue are searches for bimodal ("two-peaked") distributions and 
lattice scaling dependence on the resulting energy gap [HI [9l [lOl |TT] or scaling effects within specific heats 
(susceptibilities) P [101 ttH ttSJ [13] . For example, Ref. [10] concludes that for A ~ 1 the PT is first-order for 
[3 at and slightly above the terminal point of the phase line, but is unable to ascertain the existence of a 
tricritical point where the order of the PT would change; it would be called a critical line in (/3, k. A) space. 
Since the PT decreases in strength with increasing A, the \ = oo case presents the greatest challenge in 
the classification of the PT. Early work [8] suggested a weak first-order PT, but a recent study with large 
lattices presents evidence for a smooth crossover [13j . Furthermore, the \ = oo model augmented with 
an additional interaction leads to a line of first-order PTs which decrease in strength as the additional 
coupling approaches zero [Tl] . 

In this paper we study the phase structure of the SU (2)-Higgs model with multiple Higgs fields in the 
fundamental representation. For simplicity, we will omit Lagrangian terms which would mix more than 
one fiavour of Higgs. We find that the phase transition line extends all the way to /3 = from /? = oo for 
two or more Higgs fields regardless of the numerical value of A, in marked contrast to the case of a single 
Higgs field. We show that this enhancement of the PT is associated with overlapping PTs, as has also 
been reported for the multi-Higgs three-dimensional U{\) theory [14jFl 

The continuum S'C/(2)-Higgs action in Euclideanized space-time is 

S, = j d^x \Tr [f^'^F,,, + D^</)^t D^^^c ^ ^2^c^c ^ Ao(.^^V')'} , (1) 

where the Higgs field is in the fundamental representation of the gauge group. Including multiple Higgs 
fields in which the Higgs particles are not directly coupled to one another, i.e. they interact with each 
other only through the gauge field, gives us an action of the form. 



I n=l 



(2) 



In principle, additional quartic interaction terms are possible, including terms which couple different Higgs 
fields. In our exploratory study, we will not consider this large phase space in its entirety. Instead, we stay 
within the phase space of Eq. ([2|) . This is sufficient to observe a qualitative distinction between the single- 
Higgs and multi-Higgs theories, and is also sufficient for our discussion of overlapping PTs. The extension 
to a more complete phase space {e.g. inspired by the renormalizable scalar sector of a non-minimal standard 
model), beyond Eq. ([2]), is left for future work. 

The discretized version of Eq. ([2]) for 4-dimensional hypercube with lattice spacing a (and periodic 
boundary conditions) is 



TV 



X Kfi>iy n=l 



Kpl,n - 1)^ + Pl,n - E ^^i't'l,nUx^i(t>x+il,n) 



(3) 



For simplicity, we have chosen the couplings Ao,n so that the discretized theory contains a single four-point 
coupling A and hence the hopping parameters «:„ are defined by 

2 1 — 2A — 8k„ , , 

/^o,n = —2 • (4) 

The Higgs field has been expressed in terms of the Higgs length and an "angular component" which 
is in the fundamental representation of SU{2): 

(px = PxOLx , Px>0 , Ox e SU{2). (5) 



^By overlapping PTs we mean the influence of other Higgs fields on the PT of any one Higgs field. 
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The gauge-invariant link is defined as the gauge-Higgs coupling term with the Higgs length removed 



(6) 



which, as discussed below, is a very useful quantity for analyzing the state of the Higgs field. 

The Monte Carlo calculations we will use to investigate the nature of the Higgs phase transition employ 
a combination of heatbath steps [151 HI IIZ! and overrelaxation steps [18] to optimize the updating of the 
gauge and Higgs fields. A single Monte Carlo update for the entire system is comprised of one heatbath 
and one overrelaxation update for both the gauge and Higgs fields. The heatbath update for the Higgs 
field is implemented by generating the four real components of (j) according to the distribution 



(7) 



where 



m=l 



iTm, Tm = and Vm is from the nearest neighbour interactions. 



2 ^ 



(8) 



The distribution in Eq. ([7]) may be rearranged so that it is expressed in the form 

2 4 / \ 2l 



dp{(j)) ~ d (/i)exp < —A 



P 



1 + 



e-1 
2A 



m=l 



T 



(9) 



where ^ is an arbitrary parameter that will be determined to optimize the efficiency of the updating process. 
The four real components (j)m are each generated according to 



dp{4>m) ~ #mexp 



( 0m - -y^ 



(10) 



and the total (p is accepted with a conditional probability 



exp 



-A 



P 



1 



^-1 
2A 



The parameter is now tuned to maximize the acceptance rate jl6] : 
f + (2A - l)f - 4XC - 2A det(F) = . 



(11) 



(12) 



In cases where A and det(y) are small, approximation schemes may be used to solve Eq. (|12p . However, 
we want the Higgs updater to work reliably {i.e. with reasonable acceptance rates) for arbitrary k and 
A. Therefore, we calculate the exact positive solution for Eq. (|12p when the approximations break down. 
Also, given A > and det(y) > there is only one positive solution for Eq. (I12p . 

The phase structure of the 5^7(2)-Higgs theory can be probed by studying the way observables {e.g. 
the Higgs length, gauge-invariant link, and average plaquette) change with respect to the parameters 
K, and A. Sudden changes in the behaviour (and even discontinuities in the case of first order phase 
transitions) of these quantities characterize the location and the strength of the phase transition. While a 
variety of observables could be used to diagnose the Higgs phase transition (dynamical terms in the action 
are frequently chosen), the gauge-invariant link is our preferred quantity for two main reasons. First, it 
has a bounded range, between —1 and 1, which sets a natural scale for comparing the strength of the 
phase transition across the entire parameter space of the theory. Second, it is suitable in the extreme cases 
where A — >■ oo (where the Higgs length goes to one), and (5 ^ co (where the gauge fields go to unity and 
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Figure 1: Fast thermal cycles of the gauge-invariant link for N = 1,2,3,10 Higgs fields on a 12^ lattice 
for A = 1 in the case where all Higgs fields have identical interactions and mass terms. The solid curves 
represent (3 = 0, dashed curves represent (3 = 2, and dashed-dotted curves correspond to /3 = 4. For the 
single-Higgs case with (3 = one can see that the hysteresis effect is small compared to the statistical 
fluctuations. 

the average plaquette becomes fixed). Thus in contrast to the Higgs length and average plaquette, the 
gauge-invariant link provides useful information on the Higgs phase transition for any (3 and A. 

A typical signature of a phase transition is a hysteresis curve resulting from a heating/cooling (thermal) 
cycle. Our study does not involve temperature, but we will use the phrase "thermal cycle" to mean the 
systematic increase, then decrease, of k as follows. The gauge-Higgs coupling constant k can very loosely 
be thought of as an inverse temperature shared by the Higgs and gauge fields. By iteratively increasing 
K after a number of Monte Carlo updates, where the field configurations for the previous k are used for 
the start of the new k, and then iteratively decreasing k, a hysteresis curve in an observable (such as the 
gauge-invariant link) indicates the presence of a phase transition. The hysteresis curve can then be used 
to find the approximate location kt of the phase transition. This approach has the advantage that it does 
not require high statistics to extract information about the phase transition because it exploits the fact 
that the system does not thermalize at a phase transition with a small number of updates. This allows 
efficient exploration of a large region of parameter space for the SU (2)-Higgs model with a single Higgs 
[3l H] and is also well-suited to the multiple-Higgs case. 

For a given (/? , A), the location of the PT is initially located by performing a "fast" thermal cycle for 
the gauge-invariant link, increasing k from to 1 and then decreasing back to in steps of 0.01, using one 
heatbath and one overrelaxation update for the gauge and Higgs fields at each value of kIo A hysteresis curve 
resulting from the thermal cycle gives the approximate location of the phase transition (ht^o'^j,) fast where 
{(Tkj,) fast represents the region where the difference in the gauge-invariant link between the heating/cooling 
runs is greater than the estimated statistical fluctuations. The value of kt is then reflned by narrowing in 
on the region around the PT, k € [{kt — 2afij,)fast, (^T + 2cJkj, )/ast] i and performing a slow heating cycle 
with the same number of updates and steps in k as before, but with correspondingly smaller intervals. The 
refined location of the PT and its corresponding uncertainty is then {kt ± (^kt)sIow We emphasize that 
the signal of a PT is a hysteresis gap that exceeds the estimated statistical fluctuations within the thermal 
cycles. 



^Approximately 10 updates are used to thermalize the lattice at the initial minimum value of k and re-thermalize it at the 
maximum value. 
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Figure 2: Phase structure for the single-Higgs model on a 16^ lattice for selected values of A. The parameter 
space above the phase boundary will be referred to as the Higgs region. Note the termination of phase 
lines at small (3 for sufficiently large A, implying an analytic connection between the Higgs and confinement 
phases. 

To validate our methodology we first consider the single-Higgs case. The results for (fast) thermal cycles 
of the gauge-invariant link on a 12^ lattice for A = 1 and selected values of (3 are shown in Figure [H Figure [2] 
presents our results for the phase structure for the single-Higgs case on a 16'* lattice. The termination of 
the phase line for A = 1 is associated with the absence of a hysteresis curve for /? = and the presence of 
one for /? = 2 in Fig. [TJ For later reference, we will refer to the region of parameter space above the phase 
boundary as the Higgs region. Figure [2] compares favourably with the Ref. [1] results on smaller lattices, 
and the analytic connection that exists between the phases is expected from Ref. [2]. In addition, we have 
reproduced (not shown) the bimodal distribution functions of Ref. [10] for the gauge-invariant link used as 
a criterion for studying PTs P [9l [TOl [TT] . 

We now consider the symmetric multi-Higgs case of Higgs fields where k„ = k for all n (i.e. the 
discretized theory has a global permutation symmetry, Sn)- As illustrated in Figure [H when N > 2 
hysteresis effects exist even for (/3 = 0, A = l)o This implies a completion of the phase boundaries for the 
multi-Higgs case as shown in Figure El in contrast to the single-Higgs case where the Higgs and confinement 
regions maintain an analytic connection. The volume dependence of this conclusion has been explored by 
comparing the phase structure for 12^, 16^ and 32^ lattices; there would be no discernable difference 
between these lattices within Figure [3j 

The dependence of the PT strength on N can be assessed qualitatively from the (/3 = 0, A = oo) 
thermal cycles on a 16^ lattice shown in Fig. [H The hysteresis curve commences for N = 2 and becomes 
more pronounced as increases, indicative of increasing strength of the PT with increasing Figure [5] 
demonstrates that the signal of a first-order PT, which is difficult to ascertain in the single-Higgs case 
[7} [8| [9l [TOl II 1^ 1121 113j . becomes increasingly evident with increasing in the A = oo metastability curves 
for the gauge-invariant link obtained from hot/cold starts of the system at the transition point {k = kt)- 

Thus the PT of the symmetric multi-Higgs SU (2)-Higgs model is enhanced for small /3 as additional 
Higgs fields are incorporated into the model. To explore the possibility that this enhancement originates 
from overlapping phase transitions as found in the three-dimensional U{1) model [H], we consider an 

^For A'' Higgs fields in the symmetric case, we calculate the average of the gauge-invariant link over all fields. We have 
verified that the expectation value of the gauge-invariant link is identical for all fields as expected from the Sn permutation 
symmetry of the discretized theory. 

*The same conclusion can be extracted from Fig. [T] which also shows that the PT strength is largely unaffected for larger 

/3- 
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Figure 3: Phase structure for the symmetric double-Higgs, triple-Higgs, and ten-Higgs models on a 16^ 
lattice for selected values of A. Note the continuation of phase lines down to /3 = in contrast to the 
single-Higgs model. 

asymmetric two-Higgs case where k^. Figure E] shows the phase structure extracted from the gauge- 
invariant link for each Higgs field in the case of a large asymmetry ki = 2k\. There is an alignment of 
the PTs for A = 1 and A = oo in the /3 region where the phase lines terminate in the single-Higgs model. 
Furthermore, for A = 1 and A = cxd in the non-aligned region the lower phase line (corresponding to the 
Higgs field with larger k) closely resembles the single-Higgs case, while the upper phase line (corresponding 
to the Higgs field with smaller k) is similar to the large /3 single-Higgs model. This behaviour suggests 
that the first field to enter the Higgs region does so similar to a single-Higgs system, and influences the 
FT of the second field. Figure [6] also shows the phase structure extracted from the gauge-invariant link for 
each Higgs field in the case of a smaller asymmetry ki = I.Iki. At this level of asymmetry the phase lines 
are much more closely aligned, and the phase structure is virtually indistinguishable from the symmetric 
two-Higgs case. It should be noted that the full theory would experience multiple PTs in the non-aligned 
regions of parameter space. 

The strength of the PTs in the asymmetric two-Higgs model can be examined through the A = 1 
hysteresis curves in Figure [7| in comparison with Figures [1] and [H In the small asymmetry case, it is 
evident that the PTs for each Higgs field are comparable in strength to the symmetric two-Higgs case. 
For a large asymmetry, the first (larger k) Higgs field to enter the Higgs region resembles the single-Higgs 
system with a slight enhancement of strength, although this enhancement still results in a weaker PT than 
the symmetric two-Higgs case. However, the PT of the second Higgs field entering the Higgs region in the 
large asymmetry case behaves like a cold (large /3) single Higgs systemJl This point is further illustrated in 
Figure El where very slow thermal cycles for the gauge-invariant link are performed for A = 0.1 for the two 
different Higgs in the large asymmetry case, the two-Higgs symmetric case and the single Higgs system. 
The very slow thermal cycles were performed with 10 updates (one heatbath and one overrelaxation) to 
thermalize the lattice, 20 updates which are used to calculate expectation values of the gauge-invariant 
link, and (small) steps of 0.001 in k. In particular, the /3 dependence of the gauge-invariant link is greatly 
diminished for the second Higgs entering the Higgs region, and for (5 = 2 the first-order PT signal is absent 
for the second Higgs. This behaviour can be understood as a cooling of the gauge fields as a result of 

^In making this comparison it is important to take into account the different scales for k in Figs. [T] and [7] 
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Figure 4: Fast thermal cycles of the gauge-invariant link for N = 1,2,3, 10 Higgs fields on a 16^ lattice 
for (/3 = 0, A = oo), in the case where all Higgs fields have identical interactions and mass terms. For 
the single-Higgs case with /3 = one can see that the hysteresis effect is small compared to the statistical 
fluctuations. 
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Figure 5: Metastability of the gauge-invariant link at selected transition points for A = oo on a 16^ lattice. 



the PT of the first Higgs entering the Higgs region. Based on the comparative behaviour of the large- 
and small-asymmetry two-Higgs systems, we thus conclude that as the two-Higgs system approaches the 
symmetric limit, there will exist a mutual reinforcement of the PT strength. 

In summary, we have studied the effect of multiple Higgs fields on the phase structure of the SU (2)- 
Higgs model. The striking qualitative observation is that the analytic connection between the Higgs and 
confinement regions in the single-Higgs model is not present in the N > 2 Higgs models used in this 
work. Since our study is restricted to the action of Eq. ([3]), where all individual terms containing more 
than one Higgs fiavour are absent, it is possible that the addition of those extra terms will reveal an 
analytic connection between the Higgs and confinement regions elsewhere in that larger parameter space. 
That future exploration will allow for an interesting discussion of symmetry breaking. Furthermore, we 
see evidence for a progressive increase in the strength of the PT for increasing numbers of Higgs fields, 
and progressively stronger signals of first-order PTs as increases. Two-Higgs models with asymmetric 
hopping parameters reveal that this behaviour can be attributed to a reinforcement from overlapping PTs 
as found in three-dimensional multi-Higgs U{1) models Remarkably, a moderate 10% asymmetry 
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Figure 6: Phase structure for asymmetric two-Higgs models on a 12^ lattice for selected values of A. Error 
bars associated with determination of the transition points have been suppressed for clarity. The lower 
curve in each pair represents the phase line extracted from the gauge-invariant link for the Higgs field with 
the larger k; the upper phase line is associated with the Higgs field with smaller k. 

in the hopping parameters does not prevent this enhancement, indicating that fine-tuning of multi-Higgs 
models is not needed to realize the enhancement effects of overlapping phase transitions. 
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Figure 7: Fast thermal cycles of the gauge-invariant link for each Higgs field in asymmetric two-Higgs 
models for A = 1 and selected values of /5 on a 12^ lattice. Solid curves correspond to /3 = 0, dashed curves 
correspond to /3 = 2 and dashed-dotted curves correspond to /? = 4. 




Figure 8: Slow thermal cycles on an 8^ lattice for the gauge-invariant link(s) for the large asymmetry 
(k2 = 2ki) two-Higgs model, the symmetric two-Higgs model, and the single-Higgs model for A = 0.1 and 
selected values of /3. 
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